In this simulation study of pathogens spreading on contact networks, we sex-assortativity (r=0 to 0.6 by 0.1) and sex-trait ratios (\(\alpha\)) to see how the ratio of male-female cases changes. All simulation-parameter combinations show results from 250 replicates. Note: Sah network simulations finished for r=0, 0.6 by 0.1 but rewired networks currently running. Rewired results currently showing r=0, 0.6 by 0.2. Will update results when simulations finished.
Assorted networks were generated in two ways: with methods in Sah et al. (2019) and with a re-wiring algorithm which selectively re-wires same-sex edges until the desired level of assortativity is reached (within a small value, \(\epsilon=0.05\)). All networks had mean degree of 10.
We investigated the following sex-traits: susceptibility (SUS), infectious period (IP), and transmissibility (TRA). We investigated different strengths of these sex-traits (\(\alpha\)).
Rates of susceptibility for nodes (SUS) was modeled as
| Source -> Target | Overall transmission rate, \(\beta\) |
|---|---|
| F->F | \(\frac{2 \tau } {(\alpha + 1)}\) |
| M->M | \(\frac{2 \tau \alpha} {(\alpha + 1)}\) |
| M->F | \(\frac{2\tau } {(\alpha + 1)}\) |
| F->M | \(\frac{2 \tau \alpha} {(\alpha + 1)}\) |
Rates of transmissibility (TRA) for nodes was modeled as
| Source -> Target | Overall transmission rate, \(\beta\) |
|---|---|
| F->F | \(\frac{2 \tau } {(\alpha + 1)}\) |
| M->M | \(\frac{2 \tau \alpha} {(\alpha + 1)}\) |
| M->F | \(\frac{ 2 \tau \alpha} {(\alpha + 1)}\) |
| F->M | \(\frac{2 \tau} {(\alpha + 1)}\) |
Infectious period (INF) for nodes modeled as
| Source | Overall transmission rate, \(\beta\) |
|---|---|
| F | \(\frac{\gamma (\alpha + 1) } {2}\) |
| M | \(\frac{\gamma (\alpha + 1) } {2\alpha}\) |
Each sex-trait was modeled separately.
Sensitivity parameters tested:
Response variables:
The sex assortativity coefficient is not reported in many studies. In the Kampala network study, we found sex assortativity coefficient to be 0.26. Horton et al. compared mixing by sex in 20 studies. They measured within-sex contact for adults and children. For adults, the IQR for men was .54 to .58 and for women was .57 to .63. The plot below shows the range .54 to .63 to represent what we feel the overall range could be for adults. For children, girls had sex assorted contact while boys did not probably showing the primarily female role of taking care of children.
We think these results, both the sex assortative coefficient from Kampala and Horton’s meta analysis, show that realistic levels of sex-assortativity in populations is from 0.2 to 0.4 and we will focus our attention on these parameter ranges.
Changes in network statistics using the re-wiring algorithm to make networks of varying assortativity.
Below we examine the relationship between the epidemic threshold (\(R_0=1\)) as calculated in Kiss, Miller, & Simon (2017) and simulated epidemics on networks with different values of \(\tau\). Other parameters are given in figure caption. The epidemic threshold is given by:
\(\frac{\tau <K^2-K>}{\tau+\gamma <K>}>1\)
Comparison between analytically calculated \(R_0\) (Kiss, Miller, Simon 2017) and simulated epidemic size on non-assorted and assorted networks. Horizontal grey line shows where \(R_0=1\) and where ending epidemic size > 0. Vertical grey line approximately where \(R_0=1\). Results are consistent for SIR and SLIR models. No sex-trait heterogeneity is included in these simulated data. Other parameters: $\(\gamma=0.5; I_0=.01; \psi=0, 0.1; \delta=0.25, 1\cdot6\)
We see that the threshold matches up with simulated results (i.e., where \(R_0=1\)) epidemics are possible (i.e., outbreak size > \(I_0\)).
For \(\tau\) parameters used in this study, mean \(R_0\) on SAH networks for SLIR models is given below.
| tau | mean_r0 |
|---|---|
| 0.040 | 1.504810 |
| 0.075 | 2.452056 |
| 0.100 | 3.374769 |
While the problem investigated with these models is inspired by TB, the models themselves are very simple, and not parameterized to any specific dataset. A marked trait of TB is that, globally, 25% of the world’s population has latent infection. Here, we investigate the equilibrium latent infections and case ratio for different infection rates and strengths of sex traits. We think these data points give an idea of the parameters that we should be focusing on with these models specifically for investigating male-bias in TB. The horizontal panels are faceted by assortativity level, and in between r=0.2, 0.3 is the level of sex-assortativity observed in a previous analysis of the urban social network studied by Dr Whalen and Dr Kiwanuka (and their teams) in Uganda.
Latent equilibrium and case-ratio on Sah networks.
Only Sah results are shown here but conclusions from simulations on rewired networks were similar.
It looks like \(\tau=0.075\) or \(\tau=0.1\) results in an equilibrium latent prevalence of 25%. In the next section, we look at results with \(\tau=0.075\) and then we consider the results with a slower spreading pathogen \(\tau=0.04\) and faster spreading pathogen \(\tau=0.1\).
First, we show results for simulated male-bias with \(\tau=.075\) on SAH networks with varying \(\alpha\).
Results for the Sah networks above are representative of male-bias/assortativity results from other networks:
Results are similar on re-wired networks (only SLIR results shown below).
The graph above, which compares simulations on different network types, shows how higher male susceptibility is the only sex trait that can lead to male bias with assortativity.
With some assortativity, male bias increases with longer male infectious periods and higher male transmissibility but doesn’t seem to reach global levels of male-bias in cases.
Here we show the sensitivity of these results to different \(\tau\) values for SLIR model.
And here we show the sensitivity of these results to different \(\tau\) values for SLIRS model.
Faster spreading pathogens typically result in lower male-bias than slower spreading pathogens.
The heat plots below show which parameter combinations lead to male:female case ratios greater than 1.25 given our best guess \(\tau\) estimate of 0.075.
The effects of sex-traits and assortativity on male-bias in simulations on SAH-networks. Results shown for low transmissability pathogen. Only ratios above 1.35 are shown to visualize which parameter combinations in which models can lead to male-bias. In WHO TB report, mean M:F case report ratio is 1.77 (+/- 1 SD 0.42, range 1.35, 2.19).
The effects of sex-traits and assortativity on male-bias in simulations on SAH-networks. Results shown for low transmissability pathogen. Only ratios above 1.35 are shown to visualize which parameter combinations in which models can lead to male-bias. In WHO TB report, mean M:F case report ratio is 1.77 (+/- 1 SD 0.42, range 1.35, 2.19).
The above plot shows possible combinations that lead to male-bias (1.7, 1.9). Any combinations without color lead to male:female case ratios < 1.1.
The effects of sex-traits and assortativity on male-bias in simulations on SAH-networks. Results shown for low transmissability pathogen. Only ratios above 1.35 are shown to visualize which parameter combinations in which models can lead to male-bias. In WHO TB report, mean M:F case report ratio is 1.77 (+/- 1 SD 0.42, range 1.35, 2.19).
What if sex-assortativity has the effect of reducing clustering, path length, and increasing degree-assortativity in the real-world? Here we use rewired small-world networks, which have much different structures compared with basic small-world networks, to see how these differences are related to male-bias.
The mean peak, duration, and final size does not vary much (heat plots not shown) with assortativity.
The mean peak, duration, and final size does not vary much (heat plots not shown) with assortativity. To see the distribution of these values and results on different networks, see below box plots which are colored by assortativity values throughout this analysis.
The simulations on Sah networks suggest that assortativity does not variables much. But the other networks, which sometimes have increasing or decreasing clustering coefficients with increasing assortativity, sometimes indicate assortativity does affect final size.
Zooming in on SIR and SIRS specifically, the below plot shows how assortativity affects final size and equilibrium size specifically.
The plot above shows how different conclusions could be reached about the effects of assortativity depending on the characteristics of the network as assorativity is increased. If clustering increases as assortativity increases (SF), then increased assortativity could look to limit outbreak size. The opposite is true if clustering decreases as assortativity increases (SW).
Focusing on outbreaks on Sah networks, below we see whether differences in sex-traits impact epidemic dynamics.
In the plot above, SIR results are shown on Sah networks. They show increasing susceptibility in males relative to females can lower the overall fraction of the population that gets infected.
Finally, in SLIRS model, how does assortativity change equilibrium latent prevalence?
Increased male susceptibility relative to female susceptibility can decrease latent equilibrium infections.